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Abstract 

The neutrino oscillations in the field of a rotating deformed mass is investigated. 
The phase shift is evaluated in the case of weak field limit, slow rotation and small 
deformation. To this aim the Hartle-Thorne metric is used, which is an approximate 
solution of the vacuum Einstein equations accurate to second order in the rotation 
parameter a/M and to first order in the mass quadrupole moment q. Implications 
on atmospheric, solar and astrophysical neutrinos are discussed. 
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1 Introduction 



In the Standard Model with minimal particle content neutrinos are massless 
left-handed fermions. The question whether neutrinos have a non-vanishing 
rest mass influences research areas from particle physics up to cosmology, but 
it remains an open issue [1] . At present all hints for neutrino masses are con- 
nected with neutrino oscillation effects, namely the solar neutrino deficit, the 
atmospheric neutrino anomaly and the evidence from the LSND experiment 
[2] . Possible extensions of the Standard Model to generate neutrino masses are 
reviewed, e.g., in Ref. [3]. 

Mass neutrino mixing and oscillation in flat spacetime were proposed by Pon- 
tecorvo [4]. Later on Mikheyev, Smirnov and Wolfenstein [5] investigated the 
effect of transformation of one neutrino flavor into another in a medium with 
varying density. There have been many experimental studies exploring the 
evidence for oscillations of both atmospheric and solar neutrinos as well as 
imposing limits on their masses and mixing angle (see, e.g., Ref. [6] and ref- 
erences therein). 

The possibility to detect CP violation effects in neutrino oscillations by fu- 
ture experiments has also been considered in recent years [7,8,9,10]. Neutrino 
oscillation experiments are expected to provide stringent bounds on many 
quantum gravity models entailing violation of Lorentz invariance, so allow- 
ing to test quantum gravity theories [11,12]. Planck scale-induced deviations 
from the standard oscillation length may be observable for ultra-high-energy 
neutrinos emitted by galactic and extragalactic sources by means of the next 
generation neutrino detectors such as IceCube and ANITA [13]. Furthermore, 
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since neutrinos can propagate freely over large distances and can therefore pile 
up minimal length effects beyond detectable thresholds, there is the possibility 
to explore the presence of a quantum-gravity-induced minimal length using 
neutrino oscillation probabilities [14]. 

The effect of gravitation on the neutrino oscillations has been extensively inves- 
tigated in the recent literature, starting from the pionering work of Stodolsky 
[15]. The correction to the phase difference of neutrino mass eigenstates due 
to the spherically symmetric gravitational field described by the Schwarzschild 
metric was calculated in various papers within the WKB approximation [16,17,18,19,20] . 
The results obtained in these papers differ from each other due to different 
methods used to perform the calculation. For instance, calculating the phase 
along the timelike geodesic line will produce a factor of 2 in the high energy 
limit, compared with the value along the null line [21]. A different method was 
proposed by Linet and Teyssandier [22], based on the world function devel- 
oped by Synge [23] and defined as half the square of the spacetimc distance 
between two generic points connected by a geodesic path. Unfortunately, the 
calculation of the world function is not a trivial task. In general, it is performed 
perturbatively unless the solution of the geodesic equations is explicitly known, 
as in the very special cases of Minkowski, Godel, de Sitter spacetimes and the 
metric of a homogeneous gravitational field [24]. The effect of spacetime ro- 
tation on neutrino oscillations has been investigated in Ref. [25], where the 
Kerr solution was considered. A mechanism to generate pulsar kicks based on 
the spin fiavor conversion of neutrinos propagating in a slowly rotating Kerr 
spacetime described by the Lense-Thirring metric has been recently proposed 
[26] . Furthermore, the neutrino geometrical optics in a gravitational field and 
in particular in a Lense-Thirring background has been investigated [27]. Fi- 



3 



nally, in Ref. [28] the generalization to the case of a Kerr-Newman spacetime 
has been discussed. 

In the present paper we calculate the phase shift in the gravitational field 
produced by a massive, slowly rotating and quasi-spherical object, described 
by the Hartle-Thorne metric. This is an approximate solution of the vacuum 
Einstein equations accurate to second order in the rotation parameter a/M 
and to first order in the mass quadrupole moment 5, generalizing the Lense- 
Thirring metric. We then discuss possible implications on atmospheric, solar 
and astrophysical neutrinos. The units G — c — h — 1 are used throughout 
the paper. 



2 Stationary cixisymmetric spacetimes and neutrino oscillation 

The line element corresponding to a general stationary axisymmetric solu- 
tion of the vacuum Einstein equations can be written in the Weyl-Lewis- 
Papapctrou [29,30,31] form as 

ds2 = -/(dt-a;d(/))2 

by using prolate spheroidal coordinates {t,x,y,(l)) with a; > 1, — 1 < y < 1; 
the quantities /, cu and 7 are functions of x and y only and o" is a constant. 
The relation to Boyer-Lindquist coordinates {t, r, 9, 4>) is given by 

r - M 

t — t , X — , y — cos 9 , (f) — (f) . (2) 
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2. 1 Geodesies 



The geodesic motion of test particles is governed by the following equations 
[32] 
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where Killing symmetries and the normalization condition gapx°'x^ — — //^ 
have been used. Here E and L are the conserved energy (associated with 
the Killing vector dt) and angular momentum (associated with the Killing 
vector d,j)) of the test particle respectively, jj, is the particle mass and a dot 
denotes differentiation with respect to the affine parameter A along the curve; 
furthermore, the notation 
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(4) 



has been introduced. For timelike geodesies, A can be identified with the proper 
time by setting = 1. Let U be the associated 4- velocity vector [U -U = — 1). 
Null geodesies are characterized instead by /i = 0. Let K be the associated 
tangent vector {K ■ K — 0). 

Let us consider the motion on the symmetry plane y — 0. li y — and 
y = initially, the third equation of Eqs. (3) ensures that the motion will be 
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confined on the symmetry plane, since the derivatives of the metric functions 
with respect to i.e., fy, ujy and 7^/, all vanish at y = 0, so that y = too. 
Eqs. (3) thus reduce to 
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where metric functions are meant to be evaluated at y = 0. 



2.2 Neutrino oscillations 



The phase associated with neutrinos of different mass eigenstate is given by 
[15] 

Pt.ik)dx'' , (6) 

if the neutrino with 4- momentum P — mkU is produced at a spacetime point 
A and detected at B. 

The standard assumptions usually apphed to evaluate the phase are the fol- 
lowing (see, e.g., Ref. [33]): a massless trajectory is assumed, which means that 
the neutrino travels along a null geodesic path; the mass eigenstates are taken 
to be the energy eigenstates, with a common energy E; the ultrarelativistic 
approximation <^ is performed throughout, so that all quantities are 
evaluated up to first order in the ratio nik/E. 

The integral is carried out over a null path, so that Eq. (6) can be also written 

as 

^k= r"p^(fe)XMA, (7) 

JXa 
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where X is a null vector tangent to the photon path. The components of P and 
K are thus obtained from Eq. (5) by setting = rrik and /i = respectively. 
In the case of equatorial motion the argument of the integral (7) depends on 
the coordinate x only, so that the integration over the affine parameter A can 

be switched over x by 

= j^^ Pt^{k)^^x , (8) 
where — dx/dX. By applying the relativistic condition <^ E we find 

~ ip-a^— ^ / =dx , (9) 

2 E ^a^^x^ -1)- p{b-uj)^ 

to first order in the expansion parameter rUk/E <^ 1, where E is the energy 

for a massless neutrino and b = L/E the impact parameter. 

Therefore, the phase shift responsible for the oscillation is given by 

..,^..-^,.4.^^/" f d.. (10) 

z JxA - 1) - - a;)2 

where 

Amlj ^ml-m) . (11) 

The question as to whether neutrino oscillations should be thought of as tak- 
ing place between states of the same energy or the same momentum is still 
open. The various controversies concerning quantum-mechanical derivations 
of the oscillation formula as well as the contradictions between the existing 
field-theoretical approaches proposed to settle them are reviewed in Ref. [34]. 
The advantage of the equal-energy prescription is that the time dependence 
completely drops from the phase difference. This is also justified by the fact 
that in none of the neutrino oscillation experiments the time was measured, 
only distances between creation and detection points, as discussed by Lipkin 
and Stodolsky [36]. They showed that in the plane wave approximation 
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neutrino oscillations are experimentally observable only as a result of inter- 
ference between neutrino states with different masses and the same energy. 
All interference effects between neutrino states having different energies are 
destroyed by the interaction between the incident neutrino and the neutrino 
detector [37]. The absence of clocks in these experiments allows to consider 
the behaviour of neutrinos as a "stationary" one, only the distance between 
the source and detector being known. Therefore, the time interval is not an 
observable and only the oscillation wave length is measured. 

When the separation between source and detector is large enough, the coher- 
ence between the different mass eigenstates is expected to be lost. However, 
for atmospheric and solar neutrinos, where the source is free to move in dis- 
tances many orders of magnitudes larger, the decoherence distance will be 
even larger [38]. This is in agreement with the result quoted in Ref. [39] that 
the coherence is lost only at astronomical distances much larger than the size 
of the solar system and that this coherence loss is relevant only for supernova 
neutrinos. 



3 Neutrino oscillations in the Hartle-Thorne metric 



The exterior field of a slowly rotating slightly deformed object is described by 
the Hartle-Thorne metric [40], whose line element can be written in Lewis- 
Papapetrou form (1) with metric functions 
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and (7 = V-ZW"^ — a^- Here Pz(y) and (5/(x) are Legendre polynomials of the 
first and second kind respectively, while the functions 
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correspond to the Schwarzschild solution [q = = a). The connection between 
this form of the metric and the original one as derived by Hartle and Thorne 
was discussed in Ref. [32]. 



When expressed in terms of the standard Boyer-Lindquist coordinates (2) the 
phase shift (10) becomes 
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where terms of the order of q{a/M) as well as higher order terms in the weak 
field expansion M/r <C 1 have been neglected. In the limiting case of vanishing 
quadrupole parameter (g = 0) Eq. (14) reproduces the results of Ref. [25] for 
the slowly rotating Kerr spacetime. It is also useful to replace the parameters 
a/M and q by the Hartle- Thorne angular momentum J' and mass quadrupole 
moment Q according to 
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For instance, for 6 = Eq. (14) gives 
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The monopole term is the dominant one, due to the large distance between 
source and detector. However, describing the background gravitational field 
simply by using the spherically symmetric Schwarzschild solution is not sat- 
isfactory in most situations. In fact, astrophysical sources are expected to be 
rotating as well endowed with shape deformations leading to effects which 
cannot be neglected in general. The modification to the phase shift induced 
by spacetime rotation has been already taken into account in Ref. [25]. We 
will estimate below the contribution due to the quadrupole moment of the 
source with respect to that due to spin for solar, atmospheric and astrophysi- 
cal neutrinos by evaluating the ratio 
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In the case of the Sun we have M.q ~ 1.5 x 10^ cm, Rq a; 7 x 10^° cm and 
Jo 5 X 10^ cm^, so that 

:^~2xl0-^ |0~_4xlO-^^ (19) 

Rq Rq Rq 
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where the mass quadrupole moment has been evaluated through the rela- 
tion Qq = -J2MqRI, with J2 2 X 10~^ (see, e.g., Ref. [41]). Therefore, 
[Jq/Mq]/M% pa 5 X 10"^ and Qq/Mq fti -4.4 x 10^, so that the quadrupole 
contribution is 10^ times greater than that due to spin. 

In the case of the Earth the quadrupole term is even more dominant. In fact, 
since ^ 0.45 cm, i?® 6.4 x 10^ cm and J© ?a 1.5 x 10^ cm^, we have 



^^7x 10-^° , ^^3.7x 10-^6 , % ~ -7 X 10-13 , 

Rff, R(n Rffi 



with Qe = -J2M(sR^ and J2 ^ 10 ^, implying that [J^/M(b]/M% ^ 5.5 x 
10^ and Q(b/MI -2 X 10^^ leading to $[,7^'^V*ii^^ ~ 5 x 10^ 

The contribution due to rotation is expected to be comparable with that due 
to the deformation in the case of rotating neutron stars. Laarakkers and Pois- 
son [42] numerically compute the mass quadrupole moment Qns for several 

equations of state. They found that for fixed gravitational mass A^nS) the 
quadrupole moment is given as a simple quadratic fit, i.e., 

yviNs 

where i7ns is the angular momentum of the star and /c is a dimensionless 
parameter which depends on the equation of state. It varies between A; ~ 2 for 
very soft equations of state and A; ~ 8 for very stiff ones, for a typical mass 
of A^NS — 1.4M0. Therefore, in this case the spin term in Eq. (16) is of the 
same order as the quadrupole term. 

Finally, it is worth to note that the monopole term in Eq. (16) is always the 
leading one. In fact, for solar neutrinos the factor Mq/vbTa turns out to be 

about 2.3 X 10^^^, assuming ta ~ t© and tb ~ R® + d ~ d, being d ^ 1.5 x 10^^ 
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cm the Earth-Sun distance. The first correction to the monopole term is thus 
of the order of 10~^. For atmospheric neutrinos M.%/rBrA ~ -^e/-^® ~ 
5 X 10~^^, implying that the first correction is of the order of 10~^. 

4 Concluding remarks 

The issue of the interaction of neutrinos with gravitational fields is timely and 
has a lot of implications in astrophysics and cosmology. Neutrinos are usually 
generated by radioactive decays or nuclear reactions such as those occurring in 
the Sun, stars, accelerators or nuclear reactors and in particular when cosmic 
rays hit atoms. Neutrinos are also expected to carry off the largest amount 
of energy of an exploding star in a supernova. Neutrinos from core collapse 
supernovae can be emitted from a rapidly accreting disk surrounding a black 
hole and have been suggested to be responsible for the coohng process [43] . In 
general the coalescence of compact objects produces hot disks. The neutrino 
fiux from these objects is so large that it would be easily detected by currently 
online neutrino detectors [44] . 

Neutrino detectors of increasing sensitivity require a better understanding of 
propagation and oscillation properties of neutrinos in the neighborhood of 
massive astrophysical objects. In particular, it has been suggested that the 
gravitational oscillation phase might have a significant effect in supernova 
explosions [45]. 

We have evaluated the correction to the phase difference of neutrino mass 
eigenstates due to the gravitational field produced by a massive, slowly rotat- 
ing and quasi-spherical object, described by the Hartle-Thorne metric, general- 
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izing previous results. This is an approximate solution of the vacuum Einstein 
equations accurate to second order in the rotation parameter and to first order 
in the mass quadrupole moment of the source, generalizing the Lense-Thirring 
metric. The large distances covered by neutrinos from the emitting source to 
the detector allow to apply the weak field limit in the calculations, so that the 
Hartle-Thorne metric is enough to account for the leading gravitational effects. 
We have shown that, apart from the monopole term which is the dominant 
one, the contribution due to the quadrupole moment is much greater than 
that due to spin in the case of solar neutrinos and even more for atmospheric 
neutrinos. The effect of rotation is instead expected to be comparable with 
that due to deformation in the case of rotating neutron stars. 



Our result can also be relevant in the context of CPT and lepton number 
violation in the neutrino sector induced by gravity (see, e.g., Ref. [46] and 
references therein). In fact, it has been recently argued that the generation of 
a neutrino asymmetry can arise in accretion disks around rotating black holes 
or more generally when neutrinos are propagating in background spacetimes 
which are not spherically symmetric [47] . The proposed CPT violation mecha- 
nism leading to neutrino oscillations would be due to the different spin-gravity 
coupling of neutrinos with respect to anti-neutrinos. Of course, this is only one 
of all possible scenarios. Observable neutrino oscillations may result from a 
combination of effects involving neutrino masses and Lorentz violation [48]. 
Similar effects can also arise from violations of the equivalence principle [49] . 
Both current and future experiments on neutrino oscillations are expected to 
clarify the interplay of gravity and neutrino physics. 
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